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^ ! ABSTRACT 

' Context. A 5-night asteroseismic observation of the F8V star HD 203608 was conducted in August 2006 with HARPS, 

followed by an analysis of the data, and a preliminary modeling of the star (Mosser et al. 2008). The stellar parameters 
, were significantly constrained, but the behavior of one of the seismic indexes (the small spacing Suoi) could not be 

fitted with the observed one, even with the best considered models. 

Aims. We study the possibility of improving the agreement between models and observations by changing the physical 
properties of the inner parts of the star (to which Svoi is sensitive). 

Methods. We show that, in spite of its low mass, it is possible to produce models of HD 203608 with a convective core. 
No such model was considered in the preliminary modeling. In practice, we obtain these models here by assuming some 
extra mixing at the edge of the early convective core. We optimize the model parameters using the Levenberg-Marquardt 
algorithm. 

Results. The agreement between the new best model with a convective core and the observations is much better than 
for the models without. All the observational parameters are fitted within 1-a observational error bars. This is the first 
observational evidence of a convective core in an old and low-mass star such as HD 203608. In standard models of 
low-mass stars, the core withdraws shortly after the ZAMS. The survival of the core until the present age of HD 203608 
provides very strong constraints on the size of the mixed zone associated to the convective core. Using overshooting as 
a proxy to model the processes of transport at the edge of the core, we find that to reproduce both global and seismic 
observations, we must have aov = 0.17 ± 0.03 H v for HD 203608. We revisit the process of the extension of the core 
lifetime due to overshooting in the particular case of HD 203608. 
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■^j- , 1. Introduction to their inertia, and generate a region of mixing beyond the 

rfS ^ • t t edge of the core. This phenomenon has been investigated by 

. For mam sequence stars massive enough to show a con- geveral authorg ( see |Zahn 199ll and references therein), but 

(N ; vective core (M > 1.1M Q for solar-like metallicity) , the nQ satisfying theore tical or numerical description have been 

^ , associated mixed region plays the role of a reservoir for nu- d In practice , this region is modeled as an adiabatic 

t-J ■ clear reactions. The evolution pace of these stars and the layer aboye the cor6; whoge thickness is a fraction am of 

time they spend on the mam sequence depend directly on ^ pregsure gcale heigM „ ^ = ^ „ ) &nd where the 

the size of this reservoir. The imprecise knowledge we have elements are mixecL while it is admit ted that a ov is only 

of the mixing processes, particularly at the boundary of the & crude accQunt for the complex proce sses of mixing at the 

core, generates large uncertainties on the extension of the boundary of the conve ctive core, it is convenient and usual 

. mixed core and subsequently on the stellar age and mass for ^ the modeli of stellar interiors to adopt a repr esenta- 

■ a given set of surface parameters. Among the processes of tion of these processes depe nding on this parameter alone, 

transport of chemical elements that could contribute to the Differ ent studies have led to a wide range of a ov : between 

creation of a mixed zone beyond the edge of the convective zerQ fl Langer gjggi and about 2 jXiong 1985| ). In fact, it 

core, overshooting is the one invoked most often. fa cur l ent ly admitted that different values of a ov might be 

In the deep interior, convective elements rise adiabati- nccded tQ model stars of diffcrcnt masses and ages (see 
cally. They are accelerated until they reach the position of | C laret 2007|) . We therefore do not have precise knowledge 
convective stability, i.e. V rad = V ad . Then, the buoyancy of the amount of mixing at the edge of the corej and it is one 
forces cause a braking of the eddies m the radiative region. of the main goalg of aste roseismology to constrain it with 
It is however unlikely that they should stop abruptly at observations ( see iLebreton et1P99 fi. Mich el et al. 2 006). 
the boundary between the two regimes. They might pene- 
trate, over a distance d ov , in the regions of stability owing For intermediate-mass and hi gh-mass stars d ov is ad- 

mitted to play an important role (jMaeder 1976ft . The case 

Send offprint requests to: S. Deheuvels of low-mass stars is not as clear. When they reach the 

e-mail: sebastien.deheuvels@obspm.fr ZAMS, these stars present a small convective core that dis- 
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appears almost immediately. It has already been mentioned 
that an extra mixing at the edge of this early convective 
core might increase its longevity, by providing first more 
12 C, and then more 3 He in the center (Roxburgh 1985). In 
the specific case of the Sun, core overshooting was added 
in the models, but it was concluded that it had no rele- 
vant impact on the Sun's present structure, unless we add 
an unreasonable amount of extra mixing. Later on, stellar 
models of low-mass stars suggested that the overshooting 
at the edge of the core could make it survive almost until 
the end of the main sequence, although the phenomenon 
was not explained (Mowlavi 1993). 

In this article, we revisit this phenomenon in the case 
of HD 203608, a low- mass F8V star that presents solar-like 
oscillations. It was observed with the high-resolution spec- 
trometer HARPS at the ESO 3.6-m telescope in August 2006 
(Mosser et al. 2008, hereafter M08). The authors analyzed 
the oscillation spectrum and identified 15 £ = and £ = 1 
eigenmodes. They found a model that agrees with the phys- 
ical parameters and all the seismic parameters but one: the 
behavior of the small spacing 5isqi — v n ,Q~ (f n ,i + i / n-i,i)/2 
with frequency. Since HD 203608 is a low-mass star (less 
than 1 Mq), the effect of core overshooting was neglected 
in the preliminary modeling performed in M08. For this 
range of mass, stars are not expected to have a convective 
core on the main sequence, except for a small one which 
disappears shortly after the ZAMS. In Sect. [2J we show 
that in the case of HD 203608, with a reasonable amount 
of mixing, the early convective core can survive until the 
present age. The agreement between observations and mod- 
els is then greatly improved. In Sect.[3j we explain why the 
convective core of HD 203608 can survive, even when the 
burning of 3 He is no longer capable of sustaining it. 



2. Modeling of HD 203608 

2.1. Results from previous modeling of HD 203608 

For their modeling of HD 203608, M08 adopted the 
following stellar parameters: T e g = 6070 ± 100 if 
(jJonsell et al. 2005 1 Idel Peloso et al. 2005|) . L = 1.39 ± 
0.13 Lq (from the Hipparcos parallax it ~ 107.98 ± 
0.19 mas, Ivan Le euwcn 2007) and an observed surface 
metallicity of [Z/X]f = —0.5 ± 0.1 dex, from detailed anal- 
ysis (see M08). They searched for an optimal model by ad- 
justing the age of the star, its mass, the initial abundance 
of helium Yj, and the initial metallicity [Z/X]i to fit the fol- 
lowing observational constraints: T c g, L, [Z/X]f, Av, Si^qi, 
and Svq2i where Av represents a mean value of the large 
spacing, and Svqi, Svq2 are the small spacings defined as 



V n -l,2- 



(1) 

(2) 



Their best model fits the observational constraints on the 
star quite well, except for Svqi, for which significant dis- 
agreements remain (see Fig. [1}. 

2.2. A new modeling including extra mixing at the core 
boundary 

As shown in Appendix^ the small spacing 8vq\ is very sen- 
sitive to the central part of the stellar interior. The present 



Table 1. Physical and seismic parameters of HD 203608 
derived from the observations, from the best model with- 
out overshooting (model A), and from the best model with 
overshooting (model B). 



Observations 


model A 


model B 


Observational constraints 


Teff (K) 


6070 ± 100 


6095 


6068 


L/L 


1.39 ± 0.13 


1.379 


1.372 


[Z/X] { (dex) 


-0.5 ±0.1 


-0.560 


-0.548 


Av (^Hz) 


120.3 ± 0.5 


120.35 


120.25 


ai (.KT 3 ) 


-5.6 ±2.8 


-2.4* 


-3.6 


bi (mHz) 


0.96 ± 0.67 


2.85* 


0.84 


a 2 (.HT 3 ) 


-3.8 ±4.2 


-3.1 


-3.2 


b 2 (a*Hz) 


6.36 ± 0.93 


6.43 


5.97 


Free parameters 


Yi 




0.25 ±0.01 


0.26 ±0.01 


[Z/X]i (dex) 




-0.42 ±0.05 


-0.50± 0.07 


M/Mq 




0.95 ± 0.11 


0.94 ± 0.09 


age (Gyr) 




6.71 ± 0.75 


6.72 ± 0.59 


aov 




aov = 0" 


0.17 ± 0.03 


x 2 




9.1 


0.8 



quantities not fitted within 1-cr error bars 
parameter fixed in the optimization 



disagreement suggests that the way we model these inner 
regions should be reconsidered. As already mentioned, the 
models computed in M08 have no convective core. We tried 
to obtain models of HD 203608 with a convective core here 
by generating extra mixing at the edge of the core, as sug- 



gested by Roxburgh (1985) 

We performed a new modeling of HD 203608 includ- 
ing mixing beyond the boundary of the convective core 
by allowing d ov > in o ur models. We computed stellar 
models using CESAM2k (|Morel 1997|1 . and we derived the 
mode frequencies from these models with the Liege oscilla- 
tion code OLOSC. IScuflaire et al. 2008|) . We used the same 
physics as in M08 for our models, apart from the treatment 
of convect ion. We preferred to use the more realistic for- 
malism of Canuto & Mazzitelli (1991) For calibration, we 
computed a solar model with this treatment of convection 
and found qcm = 0.94. As in M08, we adopted this value 
for our modeling. Overshooting is described, as explained 
before, as an extension of the adiabatic and mixed region 
associated to the convective core, over a distance d ov de- 
fined as 

d ov = min(r c , a OY H P ) 

where r c is the radius of the convective core, Hp the pres- 
sure scale height, and a ov the overshooting parameter. The 
models computed with overshooting showed that the star 
can still have a convective core at its present age, provided 
the parameter a Q v is large enough. 



2.3. Optimization 

We looked for an optimal model fitting the global param- 
eters of HD 203608: T cS , \og(L/L Q ), and [Z/X] { (given in 
Table [1]), as well as its seismic parameters. For the seismic 
constraints, we adopted the mode frequencies obtained in 
the analysi s of M08. We used a mean value of the large 
separation Av and the coefficients of a least-squares lin- 
ear regression of the small spacings 8vq\ and &vqi- The fre- 
quency range of the identified modes is indeed small enough 



S. Deheuvels et al.: Asteroseismic study of HD 203608 



3 



to linearize them in the form: 

5v i — ai(v- vo) + h 
8u 02 ^ ai(v - v ) + b 2 , 



(3) 
(4) 



with vq the frequency of the maximum of the signal (y$ ~ 
2.75 mHz). The observational values of these coefficients 
are given in Table [T] We thus obtain a set of N = 8 obser- 
vational constraints referred to as j?° bs , i = 1,N. The free 
parameters are the stellar mass, age, initial helium abun- 
dance, metallicity, and the overshooting coefficient a v (see 
Table P. 

The best models are obtained by minimizing the \ 2 
function defined as 



N 

£ 

i=l 



_obs 



..mod 



_obs 



(5) 



where P f od , i = 1,N are the values obtained in the 
models for the observationally constrained parameters. In 
contrast to M08, we used the Levenberg-Marquardt algo- 
rithm (as described in |Miglio &: Monta lban 2005) to find 
the minimum of the \ z function, instead of computing a 
grid of models. This method is an interpolation between 
the Newton-Raphson algorithm and the steepest descent 
method. The steepest descent method is first used, ensur- 
ing a rapid approach when the minimum is far. When get- 
ting closer to the minimum, the algorithm progressively 
switches to the Newton-Raphson method for a faster con- 
vergence toward the minimum. The advantage is that one 
can find a minimum with only a few iterations. This kind of 
optimization is therefore less time-consuming than comput- 
ing a grid of models with 5 free parameters and a fine grid 
mesh. The uncertainties on the parameters are obtained 
from the covariance matrix of the standard errors in the 
free parameters. We took the best model found in M08 as 
a first guess and performed two minimizations: one with- 
out overshooting, giving best model A, and one with over- 
shooting, giving best model B. Model A is, as expected, 
very close to the one found by M08. The only difference 
lies in the error bars, which are bigger here. When us- 
ing the grid-search method, the error bars are obtained 
by finding the change in each parameter, which increases 
Xmin by 1. This approach is only correct if we can ne- 
glect the co rrelation between the different parameters, as 
explained in Bevington & Robinson (2003) Several studies 
have shown that correlations exist between the parameters 
(jBrown et al. 19941 IQzel et al. 20091) , which are taken into 
account in the Levenberg-Marquardt optimization, since we 
have access to the non-diagonal terms in the covariance 
matrix. This explains why the uncertainties were underes- 
timated in M08. 

Table [T] gives the physical and seismic parameters of 
both models. For model A, parameters a\ and b\ are ob- 
tained at 1.2 (j and 2.8 a of the observed values, respec- 
tively, causing a high value of Xmin — 9.1 (see Fig. [T]). In 
contrast, for model B, a\ and b\ are fitted within 0.8 a and 
0.2 tr, respectively. This results in a significant decrease of 
Xmin — 0.8 for model B. 




2 3 4 5 6 

Frequency (mHz) 

Fig. 1. Variations in the small spacing bv§\ with frequency. 
The squares represent the observational values of 6vqi. The 
darkened area corresponds to the linear regression of 5vqi 
derived from the observations within 1-a error bars. The 
dashed line represents the best model without overshooting 
(model A), and the full line the best model with overshoot- 
ing (model B). The frequency axis has been extended to 6.5 
mHz in order to exhibit the oscillation of 5vqi for model B. 



3. Discussion of the results 

3.1. Results of the fit 

When considering an extension of the mixed zone associ- 
ated to the early convective core induced by overshooting, 
we get a model that fits all the observational constraints 
better than within 1-a of the observed values. This decrease 
in the \ 2 value in fact stems from the survival of the con- 
vective core. Indeed, model B has a convective core that 
extends over about 3% of the stellar mass. The withdrawal 
of this core generates a discontinuity in the chemical compo- 
sition gradient, hence in the sound speed gradient (see Fig. 
(2) . It has already been established that such a discontinuity 

uencies as a function 
Provost et al. (1993)] 



induces an oscillation of the mode free 
of the radial order (see |Gough 1990[ ) 

derived the expressions of mode frequencies in the case of 
a discontinuous sound speed profile near the center, in the 
asymptotic approximation. Using the second-order devel- 
opment they propose, we obtained (see Appendix [XJ the 
following expression for Svqi: 



5v m {v) 



A 



B sin 



2vr- 



V 



(0) 



where n', A, B, and f(v) are defined in Appendix [Al With 
a discontinuous sound speed profile, the small spacing Spqi 
oscillates with a period V corresponding to the ratio be- 
tween the acoustic radius of the whole star and that of the 
discontinuity: 



V = 



dr/c 



(7) 



dr/c 



where r^isc is the radius of the discontinuity. We can see in 
Fig. [T]that 5voi indeed oscillates for model B, which was 
not the case for model A. 
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Fig. 2. Sound speed profile for the best model without over- 
shooting (model A, dotted line), and the one with over- 
shooting (model B, full line). A zoom of the curves shows 
the discontinuity at the edge of the convective core for 
model B. 

When the amount of core overshooting increases, the 
acoustic radius of the discontinuity in the chemical compo- 
sition gradient increases and the period of the oscillation 
decreases (see Eq. [7]). Therefore, the variations in Svqi be- 
come steeper. This shows that the behavior of Svqi in the 
models that have a convective core can be adjusted by mod- 
ifying the extension of the core, e.g. by changing the over- 
shooting coefficient. With an appropriate value of a v, we 
can correctly fit the behavior of 5v$\ and therefore reduce 
the x 2 value. We find a OY = 0.17 ± 0.03 for HD 203608. 
This result provides new input for the studies that seek to 
establish how convective core overshooting depends on the 
stellar mass (e.g. IClaret 2007|) . 

Stars in a range of mass half way between the presence 
and the absence of a convective core, such as HD 203608, are 
interesting objects for studying core overshooting. Indeed, 
if the seismic analysis of such a star can establish the pres- 
ence of a convective core, we obtain a firm lower bound 
for the amount of mixing at the edge of the core. In this 
case, the extension of the core can be derived from seismic 
parameters, such as Svqi, and we get a precise estimate of 
d ov . And if, on the contrary, no convective core is found, 
a strong upper limit for the extent of the region of extra 
mixing can be set. 

3.2. Survival of the core 

The convective core of HD 203608 disappears at about 200 
Myr without overshooting, and survives until about 7 Gyr 
with a ov — 0.17. We explain here this huge difference of 
longevity caused by overshooting. 

3.2.1. Initial core 

A convective core exists in main sequence models when the 
energy per unit mass e(m) generated by nuclear burning in 
the center is too high to be transported radiatively. When 
the nuclear flux is larger than the maximum flux -F cr it that 
can be transported radiatively (corresponding to the radia- 
tive flux in the case of critical stability for the Schwarzschild 
criterion, i.e. V = V a( j), convective motions transport part 
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Fig. 3. Profile of the total energy production rate e(m) in 
the center of the star at different stages of its evolution (top: 
50 Myr, middle: 500 Myr, bottom: 6700 Myr), for model A 
(dashed line) and B (full line). 

of the energy. To ensure the existence of a convective core, 
we must have high values of the flux of energy F(m) at 
low m, i.e. a high value of s in the center. This condition 
is equivalent to having a steep (negative) gradient de/dm 
in the center. Indeed, the star luminosity is almost entirely 
produced in the most central parts. If it is generated with a 
gentle gradient of e, the value of e in the center will be mod- 
erate (which is the case for model A, see Fig. [3J). Only with 
a steep gradient of s can it reach higher values (see model 
B in Fig. [3]). We therefore use the gradient of e hereafter 
as an indicator of the existence of a convective core. 

In most cases, we may assume that the energy produc- 
tion is dominated by a given nuclear reaction chain. Thus, 
e(m) is proportional to pXWT v , where p is the density, X 
the hydrogen mass fraction, T the temperature, v charac- 
terizes the temperature dependence of this specific reaction 
chain and W is the mass fraction of the input element corre- 
sponding to the reaction (see e.g. IHansen fc K awalcr 1994, 
chap. 6). We can then write 

„ dine dlnp d\nX dhxW dlnT , , 

Ve = ^ = aT~^ + a\ + 71 + • ( 8 ) 

a m m a In m a In m a in m a In m 

We then see that, to have a steep Ve, the dominant nuclear 
reactions must have high temperature sensitivity. 

For massive and intermediate-mass stars, the gradient 
of e is steep due to the strong temperature dependence 
of the CNO cycle (v ~ 20, IHansen fc Kawaler 1994)) . For 
HD 203608, when reaching the ZAMS, the temperature in 
the center is of about T e ~ 13 (where T e = T/10 6 K) in our 
models. At that temperature, the dominant reactions are 
those of the ppl chain. Since the pp reaction - 1 H(p,e + ^) 2 H 
- is the slowest, its generation rate is proportional to that 
of the whole ppl chain if it has achieved equilibrium. The 
temperature sensitivity of the ppl chain at equilibrium is 
therefore the same as that of the pp reaction v pp . Based on 
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Caughlan & Fowler (1988) we estimate v vp = 4.1 at our 
temperature. This is too weak for the s tar to have a con- 
vective core. However, as pointed out by Roxburgh (1985)1 
the abundances of 3 He and 12 C are in excess compared to 
their equilibrium values when the star enters the main se- 
quence. The temperature sensitivities of the burning of 3 Hc 
to 4 He and of 12 C to 14 N outside of equilibrium are esti- 
mated to be fHe3 = 16.7 and z/ci2 = 18-7 (Caughlan & 
Fowler 1988). We therefore obtain a Ve about four times 
steeper than for the equilibrated ppl chain. This ensures 
that an initial convective core is present. 



3.2.2. Withdrawal of the core 

To achieve equilibrium, the ppl chain requires that the ra- 
tio [ 3 He]/[H] corresponds to its equilibrium value. This ratio 
can be computed by assuming that the deuterium is con- 
stantly in equilibrium, which is a good approximation since 
the destruction of D is much faster than the other reactions 
in the chain (see |Clayton 1968 1. It is inversely proportional 
to the temperature. In the core, the temperature increases 
towards the center, and the abundance of 3 He should de- 
crease with r to achieve equilibrium. Since the elements 
are mixed in the convective core, the abundance of 3 He 
is constant with r, and the reactions are kept outside of 
equilibrium. The temperature sensitivity therefore remains 
high, which favors convection. One can say that convection 
is self-sustained here. 

However, e also depends on the abundance of the re- 
actant. As the star evolves, the abundance of 3 He in the 
core decreases since it is destroyed faster than it is created. 
Therefore, £hc3 decreases, until the flux of nuclear energy 
becomes lower than F cr n and the convective core disap- 
pears. This happens at an age of about 200 Myr for model A 
(without overshooting). When the core disappears, the ppl 
chain quickly achieves equilibrium (see Fig. 2]). The tem- 
perature sensitivity becomes that of the ppl chain, which 
causes Ve to be less steep. Besides, the elements are no 
longer mixed in the center, and in Eq. [51 dlnAT/dlnm > 
and dlnW/dlnm > 0. This also contributes to producing 
a more gentle Ve. 



3.2.3. Effect of an extension of the mixed zone associated to 
the core 

The existence of a mixed zone at the boundary of the con- 
vective core modifies the abundances of elements in the 



center, and it was already suggested by Roxburgh (1985) 
that it should increase the lifetime of the core. 

The peak we observe in the profile of the 3 He abundance 
(see Fig. @| is formed when the star reaches the ZAMS. 
It can be shown that the time required to achieve equi- 
librium is inversely proportional to the temperature (see 
|Clayton 1968[ ) . To the left of the peak and outside the con- 
vective core, the ppl chain is already in equilibrium, and the 
abundance of 3 He decreases towards the center because the 
temperature increases. To the right of the peak, the reac- 
tions are not in equilibrium yet, and the abundance of 3 Hc 
decreases towards the surface since the nuclear reactions 
get less efficient because of decreasing temperature. 

When adding a mixed zone at the edge of the core, 
we can see in Fig. @] that the abundance of 3 He in the 
core increases, owing to the peak we just mentioned. 
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Fig. 4. Abundance of 3 He in the center for two ages: 50 
Myr (left) and 200 Myr (right). The dashed line stands 
for model A (without overshooting) and the solid line for 
model B (with overshooting). 



Consequently, Ch C 3 increases. The convective core is there- 
fore bigger and survives longer than in the case without 
additional mixing. If the peak described above did not ex- 
ist, the mixing would not change the abundance of 3 He and 
would have little effect on the core lifetime. No such peak 
exists for the 12 C profile, and its abundance is small be- 
cause of the low metallicity of HD 203608. The role of 12 C 
is negligible compared to that of 3 He in sustaining the core 
in the case of HD 203608. 

When adding overshooting to our models, we indeed ob- 
serve an extension of the core's lifetime. However, this ex- 
tension is quite short. For example, with a ov =0.1, the core 
disappears at an age of about 1 Gyr. In model B, the core 
is still present at an age of about 7 Gyr with a ov ~ 0.17. 
This sudden increase of the core lifetime for a v > 0.1 is 
in fact caused by the start of the ppll and ppIII reaction 
chains, and later by the CNO chain, because of the increase 
in temperature in the center as the star evolves. If the 3 He 
has kept the convective core going until the ppll chain be- 
gins to compete, the convection prevents the ppll reactions 
from achieving equilibrium, in exactly the same way as de- 
scribed above for the ppl chain. The star starts burning Li 
through the reaction 7 Li(p,a) 4 He out of equilibrium, with 
a temperature sensitivity of ^Li7 = 10.8, which sustains the 
core. This is what happens for model B. If, on the contrary, 
the convective core has already disappeared at that time, 
there is no more mixing in the center and the ppll reactions 
achieve equilibrium without triggering convection. This is 
the case for model A. 

We can see in Fig. [5] that the ppll chain makes a major 
contribution to the nuclear production rate e(m). The im- 
portance of the CNO reactions is still small. As mentioned 
above, the energy brought by the ppll chain is almost en- 
tirely due to the burning of 7 Li outside equilibrium, with 
the other reactions contributing very little to e pp ii. If the 
star were only slightly more evolved, the CNO cycle would 
take over, causing the convective core to grow. 

3.3. Influence on the evolution 

We notice that, even though the structure of the core is dif- 
ferent for our models with and without overshooting, their 
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Fig. 5. Energy production rate e(m) in the central regions 
of the optimal model B. The dark line represents the total 
energy production rate. The grey areas correspond to the 
luminosity produced by the different nuclear reaction chains 
(light grey: ppl chain, intermediate grey: ppll chain and 
dark grey: ppIII chain+CNO cycle). 



age is very comparable (see Table Q]). This suggests that 
the survival of the convective core does not significantly in- 
fluence the global nuclear energy produced by HD 203608 
during its evolution. This constitutes a striking difference 
with higher mass models, where overshooting is known to 
have a strong impact on the evolution of the star, especially 
on its age, for a given T e g and L. For high-mass stars, since 
the temperature dependence of the dominant CNO cycle is 
large, it operates in a narrow area in the center, and the 
wider convective core can act as a reservoir. For HD 203608, 
the pp chain is dominant, and its temperature dependence 
is much lower. Therefore, the reactions take place in an 
area wider than the extent of the convective core, and the 
mixing has less effect on the evolution. 

This means that it would be hopeless to try to char- 
acterize an extension of the convective core in this type of 
star by classical fundamental stellar parameters alone (T c g , 
L, Al), as was done by Claret (2007) for higher masses. 



4. Conclusion 

We present here a modeling of HD 203608 based on the 
analysis of harps data performed in M08. Our main result 
is that, on this basis, we find strong evidence that this old 
low-mass star has a convective core. Models with convec- 
tive cores enabled us to solve the disagreement with obser- 
vations that was pointed out in M08 for models without 
convective cores, bringing the x 2 function from 9.1 to 0.8. 
All the observed parameters for HD 203608 are now fitted 
within l-cr error bars. 

In the case of our modeling of HD 203608, the value ob- 
tained for the a ov parameter (0.17 ± 0.03) is strongly con- 
strained. Overshooting was here used as a proxy to model 
the complex processes of transport at the edge of a convec- 
tive core, as is usually done in the present state of stellar 
modeling. Rather than finding a unique absolute value for 
a v, the current aim is to try to observationally determine 
which values of the a ov parameter are needed to repre- 
sent stars of different masses and evolution stages. In this 



respect, the value obtained for HD 203608 constitutes a 
valuable input for low-mass objects. 

We discussed in detail how the existence of a convective 
core in such an evolved low-mass star can be explained by 
a reasonable amount of extra mixing (modeled here as core 
overshooting) inducing the survival of the early convective 
core. For low-mass stars such as HD 203608, an early con- 
vective core exists because of the burning of 12 C and 3 Hc 
outside equilibrium. An extra mixing at the edge of the core 
increases its lifet ime, by bringing m ore 3 He to the center, 
as mentioned in Roxburgh (1985) Here, we showed that, 
above a certain amount of overshooting (a v ~ 0.15), the 
burning of 3 He out of equilibrium sustains the core until 
the ppll and ppIII reactions take over. Convection prevents 
these reactions from achieving equilibrium, and the burning 
of 7 Li outside equilibrium is currently keeping the core con- 
vective. This is a transitional phase before the CNO cycle 
takes over. 

The observation of low-mass stars can play a specific 
role in the study of extra mixing at the edge of stellar 
cores. Indeed, the presence or absence of a convective core 
in these stars, which can be established by seismic indica- 
tors such as the Svoi small spacing, strongly constrains the 
amount of mixing at the edge of the core. This stresses the 
interest in observing this type of star with ground-based 
observation campaigns or with the space mission CoRoT 
dBaglin et al. 2006[). 
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Appendix A: Expression of the small spacing Suqi 
with a discontinuous sound speed profile 

By neglecting the variations in the gravitational potential 
(Cowling approximation) and usi ng the asymptotic dev el- 
opment derived in Tassoul (1980), Provost et al. (1993) ob- 
tained the expression of the oscillation frequencies in the 
case of a discontinuous sound speed profile. The disconti- 
nuity is characterized by the parameter rj: 



V 



Pi^i PoCo 
PiCi + p C 



where the subscripts i and o describe the quantities at the 
edge of the convectiv e core in the inner and outer regions. 
Provost et al. (1993) showed that under these assumptions 
and provided r\ is small enough, the first-order asymptotic 
development of the frequencies is given by the expression 



with 

n = ; 

Av = 



sin 2tt 

2 TT 
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■1/2 
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where n is the radial order of the mode and I the degree. 
The parameter e depends on the reflective properties of 
the surface, and its variations with frequency are neglected. 
When there is a discontinuity, the eigenfrequencies present 
a sinusoidal oscillatory behavior. These oscillations have an 
amplitude of r\Av /tt and a period V, which corresponds to 
the ratio of the acoustic radius of the whole star to that of 
the discontinuity (see Eq. [7]) . 

Based on this, we can derive the following expression 
for Sisqi : 
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As can be seen in Fig. [TJ the mean value of the small spac- 
ing 5vox is not zero, as could be expected from Eq. IA.3I To 
understand this, one needs to push the asymptotic devel- 
opment of the freque ncies to the second order. Still from 
Provost et al. (1993)[ the expression for the frequencies is 
given by 
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where 



a n ,£ = 2tt 



1/2 
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£(e + i)v 3 + v 4 



The Vj for j — 1,4 are defined in Provost et al. (1993) 

They are of the same order of magnitude as J Q dc/r. The 
small spacing Svqi is therefore sensitive to rapid variations 
in the sound speed, and particularly in the most central 
regions, owing to the 1/r factor in the integrand. Using the 
second-order development of the freqeuncies (Eq. IA.3I) , the 
expression of Svqi is a bit more complicated. We have 
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This can be rearranged in the form 
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With a discontinuous sound speed profile, the small 
spacing 8vq\ oscillates with the same period V as the indi- 
vidual frequencies. In the continuous case, r\ = 0, and the 
expression of 8vqi reduces to the usual second-order asymp- 
totic development, which varies like 1/V, as is the case for 
model A (see Fig.[T|). 



